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Modeling the Gulf Stream system: How far from reality? 


Yi Chao , 1 Avijit Gangopadhyay , 1 

Abstract. Analyses of a primitive equation ocean model 
simulation of the Atlantic Ocean circulation at 1/6° hori- 
zontal resolution are presented with a focus on the Gulf 
Stream region. Among many successful features of this 
simulation, this letter describes the Gulf Stream separa- 
tion from the coast of North America near Cape Hatteras, 
meandering of the Gulf Stream between Cape Hatteras 
and the Grand Banks, and the vertical structure of tem- 
perature and velocity associated with the Gulf Stream. 
These results demonstrate significant improvement in 
modeling the Gulf Stream system using basin- to global- 
scale ocean general circulation models. Possible reasons 
responsible for the realistic Gulf Stream simulation are 
discussed, contrasting the major differences between the 
present model configuration and those of previous eddy- 
resolving studies. 

Introduction 

The Gulf Stream (GS) is one of the most energetic 
current systems in the world ocean, and is a dominant 
feature of the North Atlantic circulation. Description and 
understanding of the dynamics of the GS, the associated 
meso-scale eddies, and their interaction with the large- 
scale circulation are fundamental to improving our un- 
derstanding of the North Atlantic’s general circulation 
and its role in regulating the Earth’s climate. 

Accurately simulating the GS system has been one of 
the most challenging tasks in ocean modeling (Semtner, 
1995; McWilliams, 1996). Due to the complex interplay 
of the flow field and density stratification in the dynam- 
ics of the GS, a realistic simulation must accurately 
account for both wind and thermohaline forcings. Here- 
tofore, even moderately high resolution (1/2° to 1/6°) 
basin- to global-scale primitive equation general circula- 
tion models have not been able to accurately capture the 
structure of the GS system. A common problem in these 
previous simulations has been an unrealistic separation 
of the GS from the coast of North America: the separa- 
tion latitude occurs too far to the north and is accompa- 
nied by a much too strong anticyclonic (clockwise) 
circulation. 

This letter presents the recent success in simulating 
the GS system with a basin-scale primitive equation 
model. The horizontal resolution of this model is approx- 
imately 1/6°, i.e. between 15 and 20 km in the vicinity of 
the GS. After twenty-two years of integration the model 
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simulation shows significant improvements over previ- 
ous eddy-resolving, basin- to global-scale calculations. 

In particular, the GS separation and its subsequent down- 
stream meandering are in good agreement with the 
available observations. The vertical profiles of tempera- 
ture and velocity fields associated with the GS resemble 
those derived from ship-based measurements. 

Model and Experiment Design 

The experiment described here makes use of the Par- 
allel Ocean Program (POP) developed at Los Alamos 
National Laboratory (Dukowicz and Smith, 1994). The 
model domain covers the Atlantic basin from 35°S to 
80°N and from 100°W to 20°E. The water exchange pro- 
cesses across the artificially closed open boundaries are 
parameterized by 5° wide buffer zones in which the mod- 
el temperatures and salinities are restored toward the 
seasonal climatology. The model is formulated on a 
spherical grid with horizontal resolution of approximate- 
ly 1/6° (0.1875° in longitude and 0.1843° in latitude), and 
37 vertical levels. 

The model has been integrated for a total of twenty- 
two years, consisting of three experiments that differ 
from each other in their surface forcings and initial 
conditions. The first experiment was integrated for ten 
years starting from the January temperature and salinity 
distributions derived from the Levitus (1982) climatolo- 
gy and zero current. The surface wind stress and heat flux 
were based upon the seasonal climatology (Hellerman 
and Rosenstein, 1983; Han, 1984), and the surface salin- 
ity was restored to the Levitus climatology. The second 
experiment was initialized with the solution from the end 
of the first experiment. The wind stress was based upon 
the ensemble monthly mean fields derived from the Eu- 
ropean Center for Medium-Range Weather Forecast 
(ECMWF) analysis, while the heat and salt boundary 
conditions were the same as in the first experiment. The 
second experiment was integrated for two years. The 
third experiment was initialized with the solution from 
the end of the second experiment. The surface heat flux 
was derived from the ECMWF analysis based upon 
Bamier (1995)’s formulation, while the wind and salt 
boundary conditions were the same as the second 
experiment. The third experiment was integrated for an- 
other ten years. Results from only the third experiment 
are presented here. 
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Simulation Versus Observations 

The major success of the present study is the realistic 
GS separation from the coast of North America near 
Cape Hatteras. This is illustrated in Figure 1 of the sea 
surface height and surface current averaged over the last 
ten years of simulation during the third experiment. In 

3155 



CHAO ET AL.: MODELING THE GULF STREAM 


Sea Level (Annual-Mean) 



-80 -75 -70 -65 -60 -55 

Surface Current (Annual-Mean) 



Longitude (W) 


Figure 1. Climatological annual-mean maps of (top) sea 
surface height and (bottom) surface current simulated by 
the present model averaged over the last ten years of the 
third experiment. The contour interval is 10 cm. 

comparison with the previous eddy-resolving calcula- 
tions, the Community Modeling Effort (CME) for exam- 
ple, the present model shows a number of improvements 
in the GS simulation. The permanent, closed, anticyclon- 
ic circulation in the Mid-Atlantic Bight as seen in the 
CME simulations (see Figure 5 in Beckmann et al., 1994, 
and Figure 2 in Bryan et al., 1995) is absent in the present 
study. As one can see from the instantaneous map of sur- 
face current shown in Figure 2, the angle of the GS 
separation from the coast of North America near 35°N 
and 75°W is quite realistic. However, an excessive num- 
ber of cold core rings are found south of the GS. These 
cold core rings propagate westward, and eventually in- 
teract with the GS near the separation point. During such 
interactions, the GS is pushed to the north of its mean 
position, leading to the slight overshoot as seen in Figure 
1. Another major improvement of the present simulation 
is that the simulated GS can be identified as a well- 
defined meandering jet between Cape Hatteras and the 
Grand Banks. 

Figure 3 compares the observed and simulated GS 
path statistics. The observed GS position is subjectively 
determined from two-day satellite images of sea surface 
temperature during 1983-1986 (Lee and Cornillon, 
1996). The simulated GS position is objectively deter- 
mined based on the maximum gradient of sea surface 
height every three days during the last ten years of the 


third experiment. From Figure 3, it is seen that the simu- 
lated GS is generally within the observed envelope. The 
mean GS path is remarkably reproduced by the ocean 
model: the root-mean-square difference between the ob- 
served and simulated GS path is 0.44° (on the order of 50 
km) computed over 55°W-80°W. The biggest discrepan- 
cy for the mean GS path occurs at the trough of the first 
meander near 68°W. The observations show a permanent 
trough near this longitude. However, the simulated 
trough appears much deeper than the observations. The 
standard deviation of the GS path is relatively small near 
Cape Hatteras (35°N and 75^), 0.33° for observations 
and 0.52° for the model simulation, suggesting a very 
stable GS separation there. The standard deviation in- 
creases dramatically after the separation point because of 
the meandering, and remains relatively constant 
downstream. Averaged over the region 55°W-8(rW, the 
simulated GS path standard deviation is larger than the 
observed: 0.83” for observations and 1.14° for the model 
simulation. This large standard deviation in the simulat- 
ed GS path can be partly attributed to our path detection 
algorithms which sometimes cannot distinguish 
cold/warm core rings from the GS. 

Figure 4 compares the simulated vertical structure of 
temperature and velocity with the in situ profile measure- 
ments (Halkin and Rossby, 1985). Both the data and the 
model simulation are analyzed in the "stream coordinate 
reference frame. In the model analysis, the origin of the 
stream coordinate is co-located with the maximum sur- 
face velocity and the downstream direction is determined 
as the direction of the maximum transport between 0 and 
2 km along the 73°W transect. The observation is the av- 
erage of 16 transections near 73°W during 1980-1983, 
while the model simulation is the average of 1 6 snapshots 
from the last ten years simulation during the third 
experiment. From Figure 4, it is seen that the simulated 
temperature and velocity structure of the GS agree quite 
well with observations. The offshore tilt of the thermo- 
cline and the maximum velocity core is realistic. The 
simulated GS core speed is 125 cm s , slightly less than 
the observed one of 150 cm s . From the temperature 
section, it is clear that the model does not properly sim- 
ulate the distribution of 1 8°C mode water at 73°W . 
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Figure 2. Instantaneous map of surface current simu- 
lated by the present model on October 18 of year 16. 
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Gulf Stream Path Statistics 



Figure 3. The observed and simulated Gulf Stream path 
statistics. The black solid line represent the observed 
mean path, and the black dot lines represent the observed 
envelope. The red solid line represents the mean path de- 
rived from the present model. The corresponding stan- 
dard deviations of the Gulf Stream path are also shown. 


Discussions 

Determining the reasons responsible for the improve- 
ments observed in the present simulation compared to 
previous experiments will require a more systematic ex- 
ploration of model sensitivity and a more in-depth anal- 
ysis of the dynamics of the GS region than are possible 
here. At this point we can only outline the major differ- 
ences between the present study and previous experi- 
ments and conjecture as to their importance. 

The horizontal and vertical resolutions used in the 
present study are the highest among the existing z- 
coordinate, eddy-resolving calculations. The initializa- 
tion procedure is also unique in the present study. All the 
previous eddy-resolving calculations have been initial- 
ized with the solutions from lower resolution models that 
already have unrealistic representations of the GS near 
the separation point. In contrast, the present calculation 
was initialized with the observed climatological hydro- 
graphic structure, and spunup at 1/6° resolution. It ap- 
pears that a horizontal resolution near 1/6° may be barely 
adequate to properly resolve the GS but unable to correct 
for strongly biased initial conditions within the decadal 
time scale of typical eddy-resolving model integrations. 
The time required for a higher resolution model to correct 
for the biased initial conditions is presumably shorter 
(Bleck et al., 1995). The apparent sensitivity of the sim- 
ulation to initialization indicates that the higher resolu- 
tion and more accurate descriptions of the hydrography 
of the GS region (e.g., Lozier et al., 1994) should be used 
in future eddy-resolving calculations. 

The selection of surface forcings also affects the sep- 
aration characteristics. Certain characteristics pertaining 
to the cyclonic wind curl distribution in the Western 
North Atlantic can be related to the separation of the GS, 
and are found in the ECMWF operational analysis but 


not in the long-term climatology. A detailed exposition 
of the effect of different wind products on the simulated 
GS separation will be presented in a separate study 
(Gangopadhyay and Chao, Dynamical impact of cyclon- 
ic wind stress curl on GS separation, submitted to Dyn. 
Atmos. Oceans, 1996). 

Another new feature in the present study is the pres- 
ence of the Greenland-Iceland-Norwegian Seas and 
prognostic calculation of the exchanges between them 
and the Atlantic Ocean. Many previous calculations have 
northern boundaries near the sill in Denmark Strait. The 
inflow of deep water is crudely parameterized in a buffer 
zone through a relaxation of the solution toward the ob- 
served climatological conditions. It has been shown that 
the solution can be sensitive to the prescription of the 
hydrographic properties in this buffer zone (Doscher et 
al., 1994). It remains to be studied to what extent the 
overflow through the Denmark Strait is realistically rep- 
resented in the present study, and what impact this has on 
the GS behavior. It is not known whether the presence of 
the Western Mediterranean Sea plays any role in simu- 
lating the GS system, though this seems less likely. The 



Figure 4. The observed (a and b) and simulated (c and d) 
temperature and downstream velocity calculated in the 
stream coordinate system. The contour interval is 1°C for 
temperature, 10 cm s' 1 for velocity less than 100 cm s' 1 
and 25 cm s' 1 otherwise. 
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free- surface formulation for the barotropic flow and the 
unsmoothed bottom topography may also improve the 
solution. 0 . 

In summary, a horizontal resolution of at least 1/6 in 
both longitude and latitude seems to be a necessary con- 
dition to obtain a realistic representation of the GS. The 
additional factors of finer vertical resolution, unbiased 
initial conditions, more realistic surface forcings and bot- 
tom topography all help to improve the accuracy of the 
solution, but may only become effective when the mini- 
mum horizontal resolution is achieved. Further experi- 
ments and more in-depth analysis of the present simula- 
tion will be necessary to determine the relative impor- 
tance of these factors in obtaining an accurate simulation 
of the GS system. 


Conclusions 

This experiment represents a significant improvement 
in the simulation of the GS system compared to the pre- 
vious eddy-resolving calculations. Among many suc- 
cessful features of the simulation, this letter describes a 
realistic GS separation from the North American coast 
near Cape Hatteras, energetic meanders of the GS be- 
tween Cape Hatteras and the Grand Banks, and realistic 
cross-sectional structure of the GS. The relatively good 
qualitative agreement of the present model simulation 
with observed GS structure has allowed us to begin to 
pialfp direct quantitative comparisons with observations. 
These quantitative measures of model skill will facilitate 
a more meaningful comparison of different experiments 
and models than has heretofore been possible. 

Our preliminary comparisons have revealed some 
model deficiencies which can be further improved. For 
example, an excessive number of cold-core rings are 
found south of the GS. The simulation of the subtropical 
mode water needs to be further examined and improved. 
Many other aspects of the basin-scale circulation remain 
to be examined in the present simulation. Ideally it is im- 
portant to conduct many sensitivity experiments to 
explore the full range of model parameters, and to inves- 
tigate the model response to various external forcings. 
Given the improvements found here in simulating the GS 
system and the continuous increase of computing power, 
we believe that basin-scale eddy-resolving ocean models 
are finally approaching reality in describing both the 
time-mean large-scale and transient meso-scale 
circulation. 
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